Hysteresis loops of the magnetoconductance in graphene devices 
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We report very low-temperature magnetoconductance AG measurements on graphene devices 
with the magnetic field H applied parallel to the carbon sheet. The AG{H) signal depends on the 
gate voltage Vg and its sign is related to the universal conductance fiuctuations. When the magnetic 
field is swept at fast rates, AG displays hysteresis loops evident for different sizes and at different 
transport regimes of the devices. We attribute this to the magnetization reversal of paramagnetic 
centers in the graphene layer, which might originate from defects in our devices. 

PACS numbers: 72.80.Vp, 73.43.Qt, 75.20.Hr, 75.50.Dd, 81.05.ue 



Among the multitude of fields of interest, carbon-based 
nano-materials are promising candidates for applications 
in spintronics owing to their low intrinsic spin-orbit ef- 
fect, as well as the low hyperfine interaction of the elec- 
tron spins with the carbon nuclei. P, In this con- 
text, the transport properties of graphene in the pres- 
ence of a magnetic field have been deeply studied show- 
ing, for instance, that application of a field perpendicu- 
lar to the graphene plane induces quantum phenomena, 
such as the quantum Hall effect [s^,^, or weak localiza- 
tion, d, Q Recently, magnetotransport measurements 
have been used to characterize the transport regime in 
etched nano-devices, but relatively few studies have 
been carried out with the field applied along the graphene 
layer, although these allow to study purely spin-related 
effects dH. 

Several theoretical works have addressed the role of 
atomic scale defects, such as adatoms and vacancies, on 
the local electronic structure of graphene and few-layers 
graphene, fin ding that they can actually carry a magnetic 
moment. [l0l - [l2| Also the possibility of long-range mag- 
netic order has be en p redicted for certain distributions of 
such defects. [H, [ij While evidence of ferromagnetism 
has been reported for bulk graphite, fisj the experimental 
results on intrinsic magnetism in graphene are still con- 
troversial [l6|,[T7| also due to difficult experimental condi- 
tions. Since the total magnetic moment of a single layer 
graphene is expected to be extremely small, it has been 
suggested that a proof of magnetism in graphene can be 
accessed via magnetotransport measurements, [l^ How- 
ever, none of the results present in the literature have ev- 
idenced intrinsic magnetism in graphene through trans- 
port measurements so far. 

Here we focus on low-field {B < IT) parallel magneto- 
conductance (magnetic field along the graphene layer). 
We show that when the field is swept at fast rates 
the magnetoconductance displays hysteresis loops that 
we relate to the presence of magnetic impurities in the 
graphene layer. 

We obtain graphene fiakes by the standard mechanical 
exfoliation method from natural graphite. Thin flakes are 
optically located with respect to pre-patterned alignment 
markers on top of p-doped silicon wafer coated with 300 
nm of oxide. The effective number of layers is checked by 
micro- Raman spectroscopy. Metal contacts (Cr/Au or 
Ti/Pt) on the graphene sheets are obtained by electron 



FIG. 1: (Color online) Source-drain conductance for varying 
back-gate voltage at 0.04K. Inset: False-color scanning elec- 
tron microscope picture of the nanoconstriction presented in 
the text. The graphene layer has been colored in light green 
(gray) to enhance the contrast; the dark regions correspond 
to the substrate after O2 plasma etching of graphene. Source 
and drain part of graphene are indicated. 



beam lithography (EEL), electron-beam or Joule evapo- 
ration and lift-off. The underlying doped silicon is con- 
tacted (from either the top or the bottom) and is used 
as a back gate. By using this technique we routinely 
fabricate graphene devices with carriers mobility ^ up 
to 5,000 cm^V~^s~^. Electrical measurements are per- 
formed by using a lock-in amplifier with an applied ac 
voltage (< 100//V at 33 Hz) in a dilution fridge oper- 
ated at a base temperature of 40 mK and equipped with 
a three-dimensional (3D) vector magnet with sweeping 
rates as fast as 0.2 T/s. 

For some of our devices, we employed a second step of 
EBL followed by low-energy oxygen plasma in order to 
reduce the channel width. In this way we produce pat- 
terned ribbons, with different aspect ratios and lateral 
sizes ranging from 500 nm down to 50 nm, and nanocon- 
strictions of different sizes ranging from 50 nm down to 
10 nm. In the following we focus on results obtained with 
a constriction of size « 60 nm (see inset of Fig. [T]), for 
which the parallel magnetoconductance signal, defined 
as AG = [G{B) - G{B = 0)]/G(S = 0), was stronger 
(patterned samples usually display larger magnetocon- 
ductance, up to 10%, with respect to the nonpatterned 
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FIG. 3; (Color online), (a) Color scale plot of the hystere- 
sis SG (defined as the difference between the curve recorded 
while upsweeping the magnetic field and the curve recorded 
while the field is ramped back) vs gate voltage, (b) Corre- 
sponding zero-field conductance for the same gate region. 5G 
is maximum in correspondence of a maximum or a minimum 
in G{Vg) and its sign is reversed between them. 



FIG. 2: (Color online) Parallel magnetoconductance of the 
devices at three different gate voltages. All measurements 
have been performed at 0.04 K sweeping the magnetic field 
at a rate of 0.05 T/s. The curves have been selected to show 
that the magnetoconductance can be positive (a) , negative (b) 
or flat (c). In the (a) and (b) cases, a hysteresis is observed. 

ones, where the signal is usually below 1%). 

It is important to stress, however, that the main ef- 
fect described in this work (hysteresis loop in the mag- 
netoconductance) is not related to a specific type of de- 
vice, since it was observed in a variety of different sam- 
ples (examples are reported the supplemental material 
[19j). These include non- patterned graphene sheets, with 
conductivity G equal to a few Go (Go = 2e^/h) and 
mobility up to a few thousands cm^V~^s~^ at 40 mK; 
nanoribbons of different sizes and aspect ratios with con- 
ductivity G between 0.1 and 1 Go and mobility (mea- 
sured outside the transport gap region, when present) 
/J, < 100 cm^V~^s~"'^; and small nanoconstrictions with 
suppressed low-temperature conductivity (G typically 
less than 0.1 Go). 

The low-temperature conductance G of the nanocon- 
striction is shown in Fig. [1] as a function of the back- 
gate voltage Vg. The curve is characterized by a flat 
low- conducting region around 10 V (suggesting the in- 
complete formation of an electronic band gap) and by 
the presence of strong oscillations of « 0.1 Go in the sig- 
nal. These oscillations, well reproducible as a function of 
V^, evidence charge-carrier phase coherence phenomena 
(quantum interference) within our device. [20| 

In Fig. [2] we report the magnetoconductance curve as 
a function of an external held applied in the plane of 
the graphene device, taken at three representative gate 
voltages Vg. We notice that the low-temperature parallel 
magnetoconductance can be either positive or negative 



on the same device, depending on the value of V^. This 
effect is different from the behavior reported for magne- 
toconductance measurements with magnetic held applied 
perpendicular to the graphene sheet, for which the sign 
is predominantly found to be positive in small etched 
devices, Q while in larger flakes is dominated by weak 
localizations effects, [sl 

Remarkably, for low flelds the AG{B) curve depends 
on the field sweep rate cycle and we observe the opening 
of a hysteresis loop. This can be characterized by the dif- 
ference SG between the two curves obtained with increas- 
ing and decreasing the field. We notice that (i) the sign of 
6G depends on the sign of the magnetoconductance AG 
and (ii) SG changes sign with the magnetic-field orienta- 
tion. In this way SG is negative for B < and positive 
when i? > in the case of positive magnetoconductance 
[Fig. Ufa)], and the opposite in the negative magneto- 
conductance case (Fig. EJb)). No hysteresis is detected 
when the magnetoconductance is flat [Fig. EJc)]. 

To show the dependence of 6G{Vg) on the back gate 
voltage in more detail, in Fig. [Sja) we plot in color code 
the intensity of SG as a function of Vg for a small interval 
where several oscillations of the conductivity are present 
[Fig. [31[b)]. Clearly, SG{Vg) is maximum when G{Vg) 
presents either a minimum or a maximum, while its sign 
with respect to the field orientation is inverted between 
them. The corresponding magnetoconductance curves 
have positive and negative sign and can be represented 
by the curves plotted in Figs, ^a.) andllfb), respectively. 
At the inflection point between a maximum and a min- 
imum in G{Vg) the magnetoconductance is flat and the 
hysteresis is not observed: this situation corresponds to 
the case depicted in Fig. [2jc). 

The hysteresis strongly depends on the magnetic-field 
sweeping rate, as evidenced in Fig. Hl^a): for fast rates, 
such as 0.1 T/s, it closes up for B ^ 1 T and its inten- 
sity is higher, while it smoothly disappears as the field is 
swept slower. At a rate of 0.02 T/s the hysteresis closes 
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FIG. 4: (Color online) (a) Hysteresis loops obtained for dif- 
ferent field sweeping rates at base temperature. The loop 
is strongly rate dependent and it disappears for slow sweep- 
ing rates, (b) Evolution of the magnetoconductance curves 
as a function of temperature sweeping the field at 0.05 T/s. 
The hysteresis loop is still present up to ~ IK, becoming un- 
detectable when the magnetoconductance signal is flat and 
dominated by noise. 



at a field B ^ 0.5 T and 5G is about four times smaller 
than what is found for a rate of 0.1 T/s. The hysteresis 
becomes practically undetectable when the field is swept 
slower than '^0.01 T/s. This sets a time scale of few 
tens of seconds for the appearance of the hysteresis loops. 
For the slowest sweep rates, the reversible, i.e., non hys- 
teretic, magnetoconductance signal is still present. Thus 
we may distinguish two distinct components of the mag- 
netoconductance: the equilibrium part of the AG signal, 
and the hysteretic part, which depends on the sweeping 
rate. We notice that AG always shows the equilibrium 
values when the field crosses zero (i.e., SG ~ 0), while 
the hysteresis shows up only when the field is raised, ei- 
ther with positive or negative direction, away from zero 
at relatively high sweeping rates. The behavior is that 
a stronger magnetoconductance signal is observed just 
after passing zero field, relaxing toward the equilibrium 
curve at higher fields [Fig. HJa) and (b)]. 

The situation changes when the temperature raises, as 
shown in Fig. IH^b). As the temperature is increased, the 
magnetoconductance becomes flat and the hysteresis dis- 
appears accordingly, becoming practically undetectable 
above 1 K. Therefore at high temperatures the situation 
is similar to what is observed for Vg fixed at an inflection 
point of G(l/g) [Fig. (He)]. 

As previously mentioned, the hysteresis loop was re- 
producibly found in many different devices. In partic- 
ular, the appearance and the main features of the hys- 
teresis loop (including temperature and field sweep rate 
dependencies) do not depend on the size of the device 



(see supplemental information [19|). For all the mea- 
sured samples, no anisotropy is detected, as long as the 
field is applied in the plane of graphene. In addition, 
we found that the behavior is quite insensitive to many 
post-processing procedures, like annealing. Indeed, we 
systematically checked that the hysteresis loops do not 
substantially change after annealing in Ar: we used tem- 
peratures from 200-400° C and times from 15 min to 2 
h. Finally it is worthwhile to stress that hysteresis loops 
were observed in different apparatus and that we care- 
fully checked our setup to exclude any possible experi- 
mental artifact. 

To interpret these results, we first argue that the obser- 
vation of a nonequilibrium phenomenon (hysteresis loop 
of magnetoconductance) at the time scale of tens of sec- 
onds can be hardly explained within the framework of 
transport models referring to solely steady state quan- 
tities. Indeed, the characteristic time scales associated 
to scattering mechanisms in graphene are usually found 
to be orders of magnitude smaller The presence of 
electron- hole puddles domains, that are known to move 
under a magnetic field, [2l| cannot play an important role 
here, since the hysteresis loop has been found in differ- 
ent transport regimes, and even at relatively high doping 
level (see Fig. SI of the supplementary information (19i]). 
where the puddles are expected to play a marginal role. 


On the other hand, the observed behavior reminds the 
processes of magnetization reversal of small magnetic mo- 
ments at low temperature, psj Based on these obser- 
vations, we model our device to consist of two physi- 
cal parts: an electrical conductor (charge carriers in the 
graphene layer) that is coupled to localized magnetic mo- 
ments, as described, for instance, in Ref. [l8|. In this 
simple model, the hysteresis loops reflect the magnetiza- 
tion reversal of the localized moments, as the conduct- 
ing graphene layer detects the magnetization behavior 
through its magnetoconductance. 

We first address the origin of the parallel magneto- 
conductance in the conducting layer and its change of 
sign with the gate voltage Vg. In the nonballistic regime, 
magnetoconductance arises from quantum interference 
effects. [Ill Results in Figs. [5] and [3] show that the mag- 
netoconductance is positive when G{Vg) has a minimum, 
negative when G{Vg) has a maximum, and it is flat when 
G{Vg) has an inflection point. In addition, in Fig. Sfb) 
it is shown that the intensity of the magnetoconductance 
decreases as the temperature is raised, which is the same 
trend that we observed for the oscillations in the con- 
ductance G(Vg) signal. Thus the in-plane magnetic field 
tends to suppress the fluctuations in the conductance, an 
effect already observed in two-dimensional electron gases 
in heterostructures of semiconductors, [2l,[2^,but this is 
the first time, to our knowledge, that it is reported for 
graphene. Since orbital effects should be excluded in our 
case, we tentatively relate this effect to the spin split of 
the density of electronic states at the Fermi level of the 
conducting layer. More details depend on the speciflc 
regime considered as, for instance, discussed in Refs.[23| 
or [l^ , but they will bring the discussion beyond of the 
scope of this work. 
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The appearance of local magnetic moments in 
eraphene has been widely studied theoretically fiol [isl . 
|28| and, more recently, also experimentally, |16l Il7| and 
it was ascribed to the presence of defects. These rnay 
include vacancies, [3 metals or H atoms, [lil, [H [H 
Magnetic moments are also expected to be present at the 
edges of graphene. [l^ Concerning our devices, since the 
observed behavior is essentially insensitive to annealing 
process, it is unli kely that such impurities are dominated 
by H adatoms. 29] We also exclude an important role 
from edge defects, since the hysteresis is observed also in 
nonetched graphene layers and changes only very little 
with the patterning process. Concerning the presence of 
metal impurities, while we can not completely exclude 
the possibility of unwanted contaminations, a series of 
experimental facts make us confident in excluding this 
hypothesis as the dominant source of the magnetic sig- 
nal. Indeed we did not observe any change when employ- 
ing different solvents or solvents from different batches. 
In addition, we prepared our samples employing two in- 
dependent procedures to clean the Si02 substrate be- 
fore the deposition of graphene, namely O2 plasma and 
Piranha cleaning, without finding variations in the sig- 
nal. Also a possible role of the electrical contacts has 
to be excluded, since, as already mentioned in the first 
part of this Rapid Communication, different metals have 
been employed without finding any change in the signal. 
Therefore we are led to conclude that structurally intrin- 
sic defects, such as vacancies, provide the dominant 
contribution to the formation of local magnetic moments. 
In Ref. [1^, Rappoport and co-workers calculated the 
magnetoresistance of a graphene in the presence of mag- 
netic impurities, considering the case of a specific type 
of magnetic moments (H adatoms) in the variable range 
hopping regime. Although the behavior we observed is 
more general, we compare their prediction on the mag- 
netoresistance to estimate the density of the magnetic 
impurities. According to their Eq. (6), we find that the 
difference between the magnetoresistance at i? = and 
high field is related to the total magnetization Ms of the 
impurities, and hence on their total number. Comparing 
their predictions of Fig. 7 (where they found a magnetore- 
sistance drop of 60%) with our results (changing from 1% 
to 3% for the highest sweep rates), we found that in our 
case the impurity concentration is therefore « 100 - 200 
times smaller, giving an approximate concentration of 
250-500 ppm. Interestingly, this rough estimation is in 
the same order of magnitude as that found in Ref. [13] • 

We now turn our attention to the origin of the mag- 
netic hysteresis. Localized magnetic moments induced by 
defects in graphene are not expected to possess magnetic 
anisotropy, due to the very small spin-orbit interaction. 
Yet, the presence of ripples can also increase spin-orbit 
coupling, giving rise to a sizable magnetic anisotropy. 
In addition, it has been theoretically shown that local- 



ized spins may interact among them through Ruderman- 
Kittel-Kasuya-Yosida interaction mediated by free elec- 
trons, jl^j and eventually give rise to ordered magnetic 
structures. 

At very low temperatures magnetic moments may lack 
phonons to relax with, maintaining their magnetization 
even when the external field is reversed. This process, 
giving rise to a typical butterfiy-shaped hysteresis loop 
and happening at the time scale of seconds, becomes visi- 
ble at a field sweep rate of a few tens of mT/s and rapidly 
disappears for longer time scale. The phenomenon is 
known as phonon bottleneck and it has been studied for 
small spin systems. For S* = 1/2 impurities, as pre- 
dicted by the majority of theories on individual vacan- 
cies, we expect the hysteresis to rapidly disappear as the 
temperature is increased or a transverse magnetic field is 
applied. However, Fig. |H[b) shows that the magnetocon- 
ductance hysteresis of our graphene devices survives at 
temperatures as high as 800 mK. In addition, we found 
that a transverse field larger than 0.5 T is necessary to 
completely suppress the hysteresis. Therefore we con- 
clude that the localized magnetic moments should have 
S > 1/2, consistent with the phonon-bottleneck mech- 
anism. This is also consistent with what was recently 
found by superconducting quantum interference device 
measurements. [l3| 

In conclusion, we observed hysteresis loops in the mag- 
netoconductance of graphene devices with magnetic field 
applied parallel to the conducting layer and at very low 
temperature (T < 1 K). The associated time scale is of 
the order of a few tens of seconds. The observed hystere- 
sis loops are evident for different sizes and at different 
transport regimes of our devices. We ascribe this behav- 
ior to the magnetization reversal of localized spin mo- 
ments, probably arising from defects in graphene. The 
magnetization reversal is detected via magnetotransport 
properties of the charge carriers within the graphene 
layer. Based on the field and temperature dependence 
of the hysteresis, we conclude that the spin of the local- 
ized moments is higher than S* = 1/2, in agreement with 
recent works. We believe that the design and the results 
of our experiments constitute the basis for the function- 
ing of more complex graphene-based spintronic devices. 
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